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I ■ Abstract 

"t^ ' Quadratic matrix equations occur in a variety of applications. In 

this paper we introduce new permutationally invariant functions of 
two solvents of the n x n quadratic matrix equation X^ — L\X — 
£o = 0, playing the role of the two elementary symmetric functions 
of the two roots of a quadratic scalar equation. Our results rely on 
the connection existing between the QME and the theory of linear 
second order difference equations with noncommutative coefficients. 

?T>^ , An application of our results to a simple physical problem is briefly 

CS ' discussed. 

o 

Q \ ric functions. 

>: 

■^ : 1 INTRODUCTION 

Matrix language belongs to all Sciences where multicomponent variables 
and noncommutativity enter scene in the description of the system under 
scrutiny. In such situations matrix methods allow indeed compact formu- 
lations and elegant resolutions of linear and sometimes nonlinear problems. 
As an example, consider the time development of a quantum system investi- 
gated in the Heisenberg picture. One is in general faced with a matrix, often 
nonlinear, systems of differential equations typically very difficult to handle 
due to the noncommutativity of the involved observables. Even when the 
problem under investigation enables the decoupling of this system achieving 



linear differential Heisenberg equations in a n-dimensional Hilbert space, the 
application of the characteristic equation method, so useful to write down 
the general integral in the scalar case, originates a n-dimensional matrix 
nonlinear algebraic equation whose resolution cannot unfortunately rely on 
general propositions. Differently from the case n = 1 (that is the scalar 
case), the fundamental theorem of algebra does not indeed hold, so that 
the roots (usually called solvents) of an algebraic matrix equation may exist 
or not and even stipulating their existence, their number cannot be simply 
related to the degree of the equation and could in particular be infinite. It 
is thus not surprising that the theory of algebraic matrix equations both for 
its wide applicability and as research subject aiming extending the beauti- 
ful chapter of the scalar algebraic equations, has recently received a great 
deal of attention in the mathematical literature (Gantmacher 1998, Gohberg 
2006, Horn 1999). Here we recall the analysis of the spectral properties of 
the associated matrix pencil via that of the matrix equation itself (Gohberg 
et al. 1978, Krein et al. 1978). Many research papers on the properties of a 
quadratic matrix equation, especially numerical approaches, have appeared 
over the last years (Bai et al. 2005, Butler et all 1985, Dennis et all 1978, 
Highham 1987, 2000, 2001, Shurbet et al 1974). 

Our paper investigates the second degree n-dimensional matrix algebraic 
equations possessing the following canonical form 

X^ - CiX - £o = (1) 

where the unknown X and the two, generally not commuting, coefficients 
C{),Ci belong to M„(C), the algebra of all complex square matrices of order 
n. We call this equation a right quadratic matrix equation (RQME). Its 
resolution is far from being trivial essentially because the simple resolutive 
formula holding for the quadratic scalar equation (QSE) is not in general 
applicable due to the noncommutative character of the problem. 
Here right means that the linear term in eq. ([I]) has the form CiX instead of 
XCi or CiX + XC'y defining a left (LQME) or a bilateral quadratic matrix 
equation (BQME) respectively. Right, left and bilateral matrix equations 
are particular cases of the Riccati algebraic matrix equation (Hore et al. 
1999) whose canonical form is 

XAX + BX + XC + D = (2) 

In the next section we will show that the resolution of this equation is always 
traceable back to that of a related RQME. 

Thus, in this paper we concentrate on the properties possessed by eq. 
([T]) and explore whether or at which extent they may be thought of as 
generalizations of well-known properties of elementary QSE. In this spirit we 
seek analogies and differences between the two quadratic equations, matrix 
and scalar, and, in connection with eq. ([I]), we formulate the following 
questions: 



Ql) To find the number of its solvents; 

Q2) To express its solvents, if any, in terms of its coefficients and vicev- 
ersa; 

Q3) To define symmetric functions of a pair of solvents. 
It is well known that all three questions may be satisfactorily coped with 
when the quadratic equation is the scalar one. We shall show that the matrix 
nature of the unknown X of eq. ([T]) together or not with the noncommu- 
tativity between the two coefficients as well as between a solvent and the 
same £i and Cq, determine properties of eq.([T]) with no counterpart among 
those possessed by the scalar equation. Our investigation will bring to light 
interesting relations between the solvents of eq. ([1]) and its coefficients £i 
and Cq interpretable as generalizations of the classical Girard-Newton and 
Waring formulas (Sansone 1952). 

In the section 3, devoted to the questions Ql) and Q2), we introduce 
the notation, give a short sketch of the existing literature and report ad-hoc 
built examples supporting some theoretical statements. Globally speaking, 
develop step by step convincing arguments enabling us to introduce, on a 
heuristic basis, our new definition of elementary functions of two solvents of 
a quadratic matrix equation. 

The main and novel results of this paper are constructed in connection 
with our reply to Q3) and are reported in section 4 where we introduce 
and study the properties of the permutationally invariant functions of two 
solvents of the quadratic matrix equation. 



2 Mapping a Riccati algebraic equation into a RQME 

The Riccati algebraic matrix equation (Hore et al. 1999) in the unknown 
Z e M„(C) 

ZAZ + c'iZ + Z£^i+£o = 0, £^i,A,A)GM„(C) (3) 

assuming the existence of A~^ and introducing the change of variable Z = 
YA~^ may be reduced to the following BQME in Y: 

r2 ' 



Y' + CiY + YC'i + £0 = (4) 



with £1 = £1 and C'l = A~^C'iA and Cq = CqA. It is worth noticing 
that the BQME Q, in its most general case, that is when both £1 7^ and 
C'l 7^ 0, is traceable back to the RQME in the canonical form ([T]) by furtherly 
changing the variable Y to X = —£'1 — Y , where now £1 = £1 — £'^ and 

£0 = £i£i — £o- The particular case when £1 = and £'1 7^ transforms 
the BQME (j3]) into a LQME which, by applying the above proposed change 



of variable, may be traced back to the RQME ([T]) with L\ = —C[ and 
Co = —Cq. 

In all the cases above considered the original matrix equation and its 
related RQME equation are simultaneously possible or impossible and when 
the two sets of solvents are not empty, a biunivocal correspondence may be 
established. 

It is worth noting, in passing, that in the symmetric case of BQME ( 
£i = C[), called the symmetrical bilateral quadratic matrix equation ( 
SBQME), the related right equation becomes X^ = Cf — Cqi whose solu- 
tions exist if the square root of Cf — Cq exists and are denoted by the matrix 



multivalued symbol \l C\ — Cq. 
The systematic possibility of transforming LQME, BQME and Riccati alge- 
braic matrix equations into a RQME justify why in this paper we concentrate 
on the canonical form ([1]) only, thereafter simply called QME. 

3 Matrix versus scalar quadratic equation 

In this section we give our reply to the questions Ql) and Q2) beginning 
a comparison between a QSE with a QME. 

3.1 On the number of solvents 

The fundamental theorem of algebra due to D'Alembert states that each 
scalar algebraic equation of degree n in the field C always admits at least one 
root in C. This theorem, together with the well-known Ruffini's reminder- 
factor theorem, implies that the number of solutions of the equation is ex- 
actly equal to n. The examples reported in Appendix A clearly demonstrate 
that the fundamental theorem of algebra does not hold for a matrix equation 
which, in turn, may be impossible, possible with finitely-many or infinitely- 
many solvents in M„(C). 

It is easy to persuade oneself that analogous examples for matrix alge- 
braic equations of higher degree may be given. Thus, we do not have at our 
disposal a general rule directly involving Cq and Ci to assess conclusions on 
the number of solvents of a matrix algebraic equation without solving it. 

It is of relevance to underline at this point that a resolutive algorithm 
of a QME does exist (Gantmacher 1998) and stems from a theorem stating 
that if X is any solution of eq. ([T]) then g{X) = 0, where g{\) = det Q{X) = 
det(A^^ — £iA — Cq) is a polynomial in A of degree 2n. The key for the 
resolution of eq.([I]) is to exploit the knowledge of the roots of g{X) asso- 
ciated to eq.([T]) and their respective multiplicities to build an appropriate 
change of the variable X into eq.([l]) which has the merit of reducing the 
QME into a set of independent matrix linear equations in a new matrix 



unknown T which must be nonsingular. From a practical point of view the 
main bottleneck of such a recipe is of course the determination of all the 
distinct roots of (7(A), even if the conditioned resolution of the independent 
linear system of ii? equations to seek the acceptable matrices T might result 
quite cumbersome. Details on such a procedure may be found for example 
in (Gantmacher 1998), where the method is presented for right matrix equa- 
tions of arbitrary degree. Since it is practically hopeless to give a reply to 
our questions Ql) and Q2) for a generic QME using the above outlined algo- 
rithm, other approaches have been searched to investigate on the existence 
and the number of solutions of a QME. 

We quote the following result due to Dennis (1978) 

Theorem 1 If g{X) has p distinct eigenvalues {Xi}^^^ with n < p <2n and 
the corresponding set of p eigenvectors {fj}f=x satisfies the Haar condition, 
i.e. every subset of n of them is linearly independent, then the different 

solvents of the QME are exactly ( \ , if p = 2n and at least ( ) , if 

n < p < 2n, all by the form 

S = Wdiag{fii)W-^ W=[wi,...,Wn] (5) 

where the eigenpairs {pi, Wi)l''^-^ are chosen from among the eigenpairs iXi,Wi)^^ 

ofg{X). 

When p = n, the distinctness condition of the eigenvalues is not needed, 
such that we obtain a sufficient condition for the existence of a solvent 

Corollary 1 If g{\) has n linearly independent eigenvectors, then QME has 
a solvent. 

The above method furnishes only diagonalizable solvents and in consequence 
it can fail to identify all solvents and can even produce no solvents when, in 
fact, solvents exist. 

Some studies on the existence of solvents which avoid the knowledge 
of the eigensystem Q{X) are known; for example, based on the contraction 
mapping principle technique it is possible to show that if Ci is a nonsingular 
matrix and 

4||£r'llll^r'A||<l (6) 

for any subordinate matrix norm, then at least two solvents exist ( Eisen- 
feld 1973) (it's worth to note the similarity with the positivity condition of 
the discriminant for the QSE case). A similar but more restrictive condi- 
tion was derived by McFarland, Lancaster and Rokne (1977) who, using the 
Kantorovich theorem on the convergence of Newton method, have derived 
several sets of sufficient conditions for the existence of a solvent. Moreover, 
Lancaster 1966, Krein and Lager 1978 have proved that the so called "over- 
damping condition" is sufficient and ensures the existence of at least two 
solvents of a QME, called the dominant and minimal solvents. 



3.2 Solvents and coefficients 

Putting £i = 0, eq. ([I]) becomes X'^ = Cq. Such an equation admits the 
resolutive formula X = ^/Cq formally similar to that valid for the QSE, the 
substantial difference being that now ^/Cq might not exist making accord- 
ingly the equation impossible. 

When £i ^ 0, the question of the existence of a resolutive formula 
becomes doubtful. In fact, it should have the form X = F{Cq,Ci), with 
F universal and matrix multivalued function of Cq and £i. Since it must 
reduce to the well known expression valid in the scalar case, it may be 
conjectured that F must contain operations like the square root of A = 
£f + ACq. Since this square root might not exist, one would expect in 
this case to have to do with an impossible equation. But, the following 
example shows that this argument is generally false. In fact, the equation 
2^/1 \ ^ / -1 1 \ ^ ,, , / -1 1/2 



X +2\ \X—[ 1 = possesses the solvent . 



even if the square root of the matrix A = £j + 4£o = ^ ( n n ) does not 



1 
^ 
exist. 

One might wonder whether for a QME some aspects characterizing the 
scalar case such as the commutativity between Cq and Ci and the existence 
of the square root of A = £f + ACq are sufficient to represent all the solvents 
of eq. ([T]) by means of a sort of "resolutive formula"like X = (£i + \/A)/2. 



The reply is negative due to the fact that when a specific realization v A of 
the many-valued symbol vA is not expressible in the form of a polynomial 
in A, it might occur that the commutator [£i, yA] = (i2ivA — yA^i) 
does not vanish. Thus, the formula X = (£i + vA)/2 cannot be used as 
a resolutive tool even under the two simplifying assumptions, except when, 
in addition, yA commutes with £i, which is generally not true. Consider 
indeed the equation X"^ — CiX + Cf/A = 0, where Co = —Cf/A commutes 
by construction with Ci and A = 0. Of course the square root of the 
matrix A = exists and has infinitely-many realizations that is the nilpotent 
matrices of index 2 defining the set M2 ■ All these matrices, except the null 
matrix, cannot of course be expressed in the form of a polynomial in A(= 0) 
and thus we expect that "the resolutive formula "X = (£1 + \/A)/2 does 
not hold even in this simple case. This means that matrix realizations \/0 
of vO not commuting with £1 might exist. 

Taking for example Ci= ( ^ M = £f and 7o = ( _^ -1 ) ^ ^2^^ 
where 

AAf ^{ivf = (^^ ^^^/a,b,ceC,a' + bc = 0} (7) 



we easily check that (£i + \/0)/2 is not a solvent of the equation X^ — 
I I X + I I I = 0. It is possible to prove that the only solvent 

of this equation is >Ci/2 and that it is obtainable from (£i + \/0)/2 solely 
choosing yO = 0. On contrary, the simple QME {X — E/2)'^ = admits 
infinitely-many solvents, namely S = E/2 + A'^2 > even if once again it has 
A = 0. Thus, for this example, the formula X = (£i + \/A)/2 is valid since 
Ci = E commutes with every realization of vA- Another intriguing example 


is provided by the equation X'^ — CiX = with £i = ( 1 | and 


Cq = 0. In this case [£i, £o] = 0, A = £f = | 1/2 | and the formula 



/O 
X = (£i-|-vA)/2 does not include the matrix I 1 | which, in turn, 

is a solvent of the equation. This last example is interesting because it shows 
that the resolutive formula may fail since it does not capture all the solvents 
of the QME. 

Summing up, we conjecture heuristically that a general resolutive for- 
mula for QME might not indeed exist. It is however of relevance to observe 




that when yA exists and [vA, -^i] = 0, for at least one realization vA of 



/A, then the QME admits the solvent X = (£i + VA)/2. 
By the way, it is interesting to note that vA exists when QME ([1]) admits 
a solvent S commuting with Ci. We may indeed write that S'^ — CiS — 
£.0 = 3^- CiS/2 - S£i/2 + Cj/A - A/4 = {S - Ci/2)^ - A/4 = ^ 
{S — £i/2)^ = A/4 which, by definition, assures that A admits square root. 

Reversing our point of view, we concentrate now on establishing whether 
given two matrices Si and S2 G M„(C), there always exists at least a QME 
admitting Si and S2 as solvents (eventually, among others solvents). We 
know that in the scalar case this important question has a unique positive 
reply, that is the knowledge of two roots of a scalar equation determines 
the coefficients of the equation (putting legitimately its director coefficient 
equal to 1). What can one say about in the case of QMEs? The reply to 
this question is more articulated than in the scalar case in the sense that all 
the possibilities may now take place: the given Si and 5*2 may be solvents of 
only one determinable QME; no QME could exist simultaneously admitting 
the two given matrices as solvents or there might exist infinitely-many QME 
satisfying such a condition (see Appendix B). 

In the following a matrix pair M, N e M„(C) such that det(M-iV) 7^ 
is called complete, accordingly with the definition introduced by ref. (Krein 
and Langer 1978) in the context of quadratic operator equations. 



We now may prove the following 

Theorem 2 Two matrices Si and S2 G M„(C) are simultaneous solvents 
of only one QME X'^ — CiX — Cq = 0, if and only if the matrix pair Siand 
S2 is complete and the coefficients L\ and Cq may he expressed as 

^Q = S2{Si — S2) Si — Si {Si — S2)~ S2 (8) 

£1 = isf - si){Si - S2r' (9) 

Proof : 

The existence of a unique QME admitting Si and 5*2 as simultaneous 
solvents of eq. ([I]) requires that 



Si — jCi<Si — jCq — 
S2 — C1S2 — Co = 

or equivalently that 

Cq = Si — CiSi = 5*2 — C2S2 



(10) 



Jl — D2 — -'-ll'-'l ~ '-'2J 

Regarding the second equation as a matrix linear equation in Ci we know 
by hypothesis that it has only one solution. This fact necessarily requires 
that rank{Si — S2) = n, that is {Si — S2) is not singular and thus necessarily 
that 

Ci = {Sl-Sl){Si-S2)-' (12) 

To determine Cq in the desired form we observe that 

-^0 = 5'i — {Si — S2){Si — S'2)^ ^i = 
Si{Si — 52)^ {Si — S2) — {Si — S2){Si — S2)^ Si 

= Si{Si - S2)-'Si - Sf{Si - S2)-'S2 (13) 

Reciprocally, if the pair Si and 5*2 is complete and Cq and £1 are given 
by eqs. ([8]) and ([9]) respectively, then it is immediate to see that Si and ^2 
are solvents of the corresponding QME. 

If {Si — S2) is not invertible, then the necessary condition Sf — S2 = 
Ci{Si — S2) may be incompatible with the existence of a QME having Si and 
5*2 as solvents or may give rise to infinitely-many possibilities for Ci and Cq. 
This means that there might even exist infinitely many QME possessing Si 
and 5*2 as simultaneous solvents. In the appendix B we give an example of 
two matrices ^i and 5*2 which cannot be simultaneous solvents of any QME 
and two matrices S'l and S2 which are simultaneous solvents of infinitely 
many QME. 

We underline that eqs. dS]) and ^ may be indirectly constructed exploiting 
lemma 2.14 of ref.Gohberg 1982. 



4 Symmetric functions associated to a QME 

4.1 Elementary symmetric functions associated to a QME 

It is interesting to point out that under the condition of vaUdity of eqs. 
dH]) and dS]), the further assumption that [^i, 52] = leads us, as expected, 
to the well known elementary relations Cq = —S1S2 and £1 = 5i + 5*2 
which in turn implies that [£o)>Ci] = 0. Even if we are sure that, being 
in addition A = {Si — 82)^, \/A exists, we do not fully recover all the 
properties holding in the scalar case since, as previously noted, there might 
exist matrix values of vA not commuting with Ci. The two expressions 
Si + 5*2 and 5*1 5*2 play an important role in the context of the theory of 
quadratic scalar equations. They indeed are the two associated elementary 
symmetric functions of the two (only and always existing) roots of the QSE, 
in terms of which every rational symmetric function of the same roots 5i 
and S2 is rationally expressible. It is immediate to persuade ourselves that 
for a QME admitting at least the two solvents Si and S2, matrix functions 
like Si + S2 and 5*1 5*2, although the simplest ones possessing permutational 
invariance under the exchange of 5*1 with S2, cannot occupy a special place 
for the QME simply because they are not invariant in the set S of all the 
possible solvents of the QME. This means that in general Si + S2 ^ S'i + S'2 
and/or S1S2 / 5*1 5*2 with 5i,52,5[,52 G S so that these functions are 
not biunivocally determined by the coefficients of the QME under scrutiny. 
Theorem 1 suggests to identify the right hand sides of eq. dH]) and Q as the 
elementary symmetric functions of two solvents to be used when dealing with 
a QME. These expressions indeed are not only permutationally invariant 
but their matrix values do not change when we substitute a complete pair 
of solvents Si, S2 with another one complete too. For this reason we are 
stimulated to introduce the following 

Definition 1 Let Si and S2 he a complete pair of solvents of a QME (Op. 
The two expressions 

S2(Si, 52) = {Si - Sl){Si - S2)-^ (14) 

and 

n2(Si, 52) = sl{Si - S2)-^S2 - sl{Si - S2)-^Si (15) 

are called elementary symmetric functions associated to a QME. 

4.2 Symmetric functions associated to a QME 

It is noteworthy to observe that for a given QME, simple permutationally 
invariant expressions like Sp^2 = S\ + S2 or Tip^2 = -Sf S'2 + S2S1 with 
p = 1,2,..., generally lacking of invariance in S whatever p is, cannot 
be in general expressed in terms of Cq and Ci. We, thus, claim that it 



is not true that every symmetric function of two solvents of a QME not 
possessing permutational invariance in S, is rationally expressible in terms 
of the elementary symmetric functions associated to a QME (fTH) and (f36]) 
only. In particular, this means that even Si^2 = 81 + 82 and ni^2 = 5*1 52 + 
S'2<S'i suffer such a difficulty and for this reason we call them irreducible 
symmetric functions. We are however going to show that when 81 and 82 
constitute a complete pair of solvents of a QME, then Sp^2 and U.p,2 may 
be represented in terms of the elementary symmetric functions Ti2{8i,82) 
and n2(S'i, S2) given by eqs. (fM|) and ([5H]) respectively as well as in terms 
of the irreducible expressions 81^2 and ni^2 respectively. We anticipate that 
such a representation of 8p^2 niay be interpreted as a generalization of the 
Waring formula, well-known in the context of QSE. 

To this end we exploit a recently published result concerning the exact 
resolution of the matrix Cauchy problem 



Yp+2 — ^oYp + CiYp^i 
Yo = 0, Yi=B 



(16) 



where the variable p runs in N and Yp, L\ and £0 belong to M„(C) (Jivulescu 
et al. 2007). We have been able to extend the successful technique envisaged 
in this reference to solve the more general Cauchy problem 



^+2 — l^oYp + L\Ypj^\ 
Y^ = A, Yi=B 



(17) 



getting the following always existing unique matrix solution 



Yp = ttpA + PpB, 



(18) 



where 



ar. 



and 



I3p 



rp-2 

t=0 



E 



(t) ^(p-2-2t) 



}Co 



rEnii 







(i)^(p-l-2i) 



} 



if p>2 



if p 
if p 



1 





if p>l 
if P = 



(19) 



(20) 



where E is the unitary matrix. The symbol {Cq C^ }, introduced in 
(Jivulescu et al. 2007), denotes the sum of all possible distinct permuta- 
tions of u factors £0 and v factors £1. 

As an immediate application of eqs. ()17p and (llSp . let firstly notice that 
if S is a solvent of eq. ([T|), then 8^ — Ci8 — Cq = and consequently 



10 



S"^"*"^ — CiS^'^^ — CqS^ = 0, where p G N. This means that the p— th power 
of any arbitrary solvent S of eq. ([T]) satisfies the fohowing Cauchy problem 



5P+2 = CoSP + CiSP+^ 
S° = E, S^ = S 



_ 17 ci _ e (21) 



which, coinciding with that expressed by eq. (fT7|) . leads to the following 
formula 

SP = l3pS + OpE, (22) 

It is worth noticing that the two operator coefficients Op and /3p ap- 
pearing in the right hand side of eq. ()22p are completely determined by 
the knowledge of the coefficients Cq and £1 only. Thus, eq. (|22]) can be 
regarded as an effective linearization of S^, for any solvent S, becoming in- 
deed a computational tool in the sense that the expressions Op and f3p, for 
a given p, may be calculated only once, whatever the solvent S is. 

In passing we observe that if S represents the transfer matrix of a "box" 
then, whatever its nature is, 5^ may be interpreted as the output after p 
identical iterations, the QME admitting 5 as solvent implicitly defining the 
transmission features of the box. In the last section we will exploit this point 
of view to deduce the transfer matrix of an A^-period potential. 

Eq. ()22p can be generalized. Given indeed r-simultaneous solvents 
Si, S2, ■ ■ ■ ,Sr of eq. ([T]), let's introduce 



Sp,r = Y.Sf (23) 



i=l 



a permutationally invariant quantity under the exchange of any i with any 
jj (^) J = 1) • • • 5 '')• It is immediate to persuade oneself that such a quantity 
satisfies the Cauchy problem 



So,r = rE, Si,r = Y. Si 



which, in view of eq. (fTH]) . admits the solution 

^P,r = ^py£,Si\^rapE (25) 

This formula is at glance deducible from eq. (|22p which in turn is recoverable 
from eq. (|25p simply putting r = 1. We emphasize that arbitrarily fixing 
p, the expression of Sp^r given by eq. (f25]) is linearly related to Si^r in the 
sense that the operator coefficients ap and /3p are independent both on the 
number r of solvents involved and on the group of r- solvents chosen in the 
set S of all solvents of eq. ([T]) . 

11 



We recall that the set S may be also infinite thus permitting us to chose 
in infinitely many ways the family of r- solvents satisfying eq. (|25p . In 
addition, like eq. (|22p . this equation may be regarded as a computational 
tool for the p-power sum of r solvents of eq. ([1]). 

We wish to notice that putting p = r = 2 into eq. (|25|) . under the 
condition that Si and S2 are a complete pair of simultaneous solvents of eq. 
dB), yields 

82,2 = S2(5i, ^2)81,2 - 2U2{Si,S2)E (26) 

which may be interpreted as generalization to the matrix case of the classical 
Girard- Newton formula for symmetric polynomials in two scalar variables 
(Sansone 1952). 

In the same spirit and hypothesis we are also able to get an extension to 
the matrix case of the well-known Waring formula for two variables. Putting 
indeed r = 2 into the identity (j25|) . leaving instead the integer p free to run 
in N, we immediately get 

Sp,2 = f3p{Si + S2) + 2apE (27) 

The peculiar feature of eq. (I27p is that the presence of the opera- 
tors ap and Pp amounts at rewriting the p-th power sum Sp^2 in terms 
of T,2{Si, 52), n2(5i, S2), the above defined elementary symmetric functions 
associated to the QME, for this reason we claim that eq. (j27p may be inter- 
preted as generalization to the matrix case of the Waring formula for two 
scalar variables (Sansone 1952) . 

Another direct consequence of eq. (|17p concerns the permutationally 
invariant quantity 

np,r= Yl (SfSj + SiSP) (28) 

l<i<j<r 

It is not difficult to see that the quantity {SfSj + S^Si), for any given i 
and j, satisfies eq. (f^ . that is 

SfSj + SPSi = Pp{S, Sj + Sj S,)+ap{S, + Sj ) (29) 

from which we immediately deduce that 

np,r = Pp'^i^r + {r- l)apSi,,. (30) 

Applying to eq. ([50]) the same arguments used to deduce eq. ([^7|) from 
eq. (!25]) immediately yields 

np,2 = apSi,2 + f3p'ni^2 (31) 

which is of relevance since it expresses ^p,2 in terms of the irreducible expres- 
sions Si^2 = (5'i + 'S'2) and ni^2 = S1S2 + S2S1 and the elementary symmetric 

12 



functions associated to the QME J^2iSi, S'2), n2(S'i, 52), considering that Op 
and (3p are functions of Cq = T,2{Si, S'2) and £1 = — n2(S'i, S'2). 

It is interesting to stress that the unavoidable presence of 81^2 in eq. ([271) 
and of ni^2 too into eq. ()3ip does not allow to write down the permuta- 
tionally invariant expressions Sp^2 and ^p,2 in terms of the noncommutative 
coefficients Cq and £1 only. This observation confirms our claim that the 
fundamental theorem regarding the symmetric functions of two scalar vari- 
ables does not hold in the matrix case. 

The elementary symmetric functions associated to the QME ([T|) may be 
easily generalized. To this end let's write eq. ()22p for both elements Si and 
S2 of a complete pair of simultaneous solvents 

Si = PpSi + ap , . 

Sf = PpS2 + ap 

It is simple to solve this system of equations in the unknowns Up and /3p 
obtaining 

ap = Sf (Si - S2)-^Si - Sf (Si - S2)-^S2 (33) 



and 



/3p = (Sf-S2P)(Si-S2)-^ (34) 



Eqs. dSD and dH]) reduce to -n2(Si,S2) and i;2(Si,S2) when p = 2. In 
addition, the right hand side of both equations are certainly invariant in the 
set S. For this reason, we are induced to introduce the following 

Definition 2 Let Si and S2 be a complete pair of simultaneous solvents of 
a QME |7p. The two expressions 

Sp(Si, S2) = (S? - S^)(Si - S2)-^ (35) 

and 

Up{Si, S2) = S[(Si - S2)-^S2 - S^(Si - S2)-^Si (36) 

are called symmetric functions associated to a QME. 

It is of relevance that Sp(Si, S2) and np(Si, S2) may be expressed in terms 
of elementary symmetric functions associated to a QME S2(Si,S2) and 
n2(Si,S2) as easily seen considering the definition of Op and (3p and in 
the view of the eqs. ([8]) and ([9]). If we restrict our definition of symmetric 
functions associated to a QME to all those expressions of two complete 
solvents satisfying the condition to be permutationally invariant into S then 
such expressions like (j35]) and ([36j) must be expressible in terms of Cq and 
£1 and then in terms of S2(Si, S2) and n2(Si, S2). 

Thus, we claim that any symmetric function associated to a QME may 
be written down in terms of elementary symmetric functions. 
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5 A simple application 

If the order of the unknown matrix X of eq. ([l|) is fixed at its minimum 
value 1 we recover the scalar case getting simultaneously rid of both the 
matrix nature of the problem and of any complication stemming from the 
noncommutativity between Cq and £i as well as between X and the two co- 
efficients. A simpler but not trivial problem, someway intermediate between 
the full matrix noncommutative eq. ([1]) and its scalar version, may be singled 
out requiring solely the commutativity between the two coefficients Cq and 
£i. We emphasize that such a restriction does not spoil eq. ([I]) of interest 
enabling us, in turn, to bring to light some peculiar properties consequent 
to its matrix nature only. We stress indeed that on the basis of the results 
reported in the previous sections of this paper, the hypothesis [i2o)>Ci] = 
does not lead to any effective recipe to build up the set of solvents of eq. 
([1]) which may be still empty or contain finitely-many or infinitely many 
elements. The concrete advantage provided by the condition [£o)>Ci] = is 
that it makes the matrix expressions of ap and fip algebraically manageable 
disclosing their connection with the Chebyshev polynomials of the second 
kind Up{x) (Murray-Spiegel 1998) defined as follows 



Up[x] 



[p/2] 
m=0 



-ly 



{p 



m] 



m\(j) — 2m)! 



(2x) 



p— 2m 



(37) 



Taking indeed into consideration the commutativity between Cq and Ci 
as well as the fact that the number of all the different terms appearing 
in the operator symbol {Cq Cf^ coincides with the binomial coefficient 

u + V 
m 



with m = min(n, u), we may rewrite for p >2 that 



ar. 



[(p-2)/2] 

E 

t=0 



P - 2 - t \ f^i p_2-2t 



t 



-'-0 ^1 



(38) 



while ao = E and ai = 0. Assuming the existence of the inverse of the 
square root of Cq opens the way to express Op as 



[(P~2)/2] 



Oir, 



t=0 



-{-CoY'^ Y: (-!)*( ^"'"' ](A(-/:o)-^)(^-^-^*) (39) 



that is 



P 

dp = -{-Copl(p_2 



:Cii-Co 



, P>2 



(40) 



Following a similar procedure and under the same hypothesis we have 
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[{p-l)/2] 

E 

[(p-l)/2] 



p-l-t \ ^4^p_i_2t 



t 



-o-^i 



(-£o)(^-^)/2"E' ■(-!)* (^ \ M (A(-/:o)-^)(^-^-2*) 



from which 



/3p = {-CoY-i Up-i 



lc^i-Co)-'^ 



, P>1 



while /3o = 0. Inserting eqs. (00]) and (f^T|) into eq. (f^^ yields 



5P 



-A) 






:p-l 



^^.<- 



■Cr 



S-{-Cof^Up 



\c. 



-CoY 



which reduces to the trivial identities S = S and E = E, when p 

p = respectively. 

If the two coefficients Cq and £i appearing in eq. 



(41) 



E, p>2 

(42) 
= 1 and 



commute, then the 



matrix S 



This amounts at claiming that S is also solvent 



right hand side of formula ()42|) provides the definition of the transfer matrix 
of a single "equivalent box " substituting all the original p ones. 

It is worth noting that the practical use of eq. ()42p is greatly facili- 
tated when the minimal polynomial of S has a degree 2. We remind that 
the monic polynomial of minimum degree p{x) such that p{S) is the zero 
matrix is called the minimal polynomial of S and that it always exists and 
is unique (Gantmacher 1998). Let us suppose that the distinct roots of 
the characteristic polynomial of S, a given solvent of eq. ([T]), are only two 
and denote them by Ai and A2. If in addition we stipulate that S is di- 
agonalizable, it is well-known that, then, its minimal polynomial coincides 
with the characteristics polynomial of the auxiliary two by two diagonal 
Ai 
A2^ 

of the QME S^ - {trS)S + {detS)E = 0, relatively simpler than eq. (^ 
since Cq and £1 are now proportional to E. As a direct consequence for- 
mula (I42p may be easily exploited to find S^, for any p. We underline that, 
when S is of order two, the applications of this formula may be made sys- 
tematically easier, simply having recourse to the Cayley-Hamilton identity 
S*^ — {trS)S + {detS)E = 0. Quite recently, for example, this approach has 
been used to express the A^- period transfer matrix M^ of a photonic crystal 

1/t r*/t* 
r/t 1/t* 

and r are the complex transmission and reflection coefficients of the unit cell 
respectively, satisfying the condition |rp -|- |tp = 1 (Benickson et al. 1996). 
It is easy to persuade ourselves that the Cayley-Hamilton theorem leads in 
such a case to the following simple QME 

M"^ -2cosPM + E = 



as a linear function of the unit cell-matrix M 



where t 
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where cos/3 = Re{l/t). 

Thus, resorting to eq. ()42p . we may soon estabhsh the transfer-matrix re- 
duction identity 

M^ =Un-i[cosI3]M-Un-2[cosI3]E 

which may be put into the following closed form (Murray-Spiegel 1998) 

j^N _ sinNIi ^.j- _ sin(Af-l)/3 j^ 
sin/3 sin/3 

This expression coincides with that differently derived in ref. (Bendickson 
et al. 1996). 

6 Conclusive remarks 

In this paper we have studied the quadratic matrix equation X'^ — CiX — 
£o = bringing to light new peculiar properties of its, clearly traceable 
back to the inherent noncommutative character of the problem. Driven by 
the knowledge of any aspect concerning the simple classical theory of the 
scalar quadratic equations, we have formulated and coped with some basic 
questions putting into evidence remarkable differences between the matrix 
and the scalar case. Exploring the meaning and the consequences of such 
differences we have succeeded in introducing the definition and the explicit 
expressions of the matrix elementary symmetric functions associated to a 
QME by which some notable identities involving solvents of the QME have 
been derived. 

Other questions greatly stimulate our attention like for example the no- 
tion of multiple solvents or how to define a factorization procedure, if any, 
of a quadratic matrix polynomial from the knowledge of its finitely-many or 
infinitely-many solvents. We believe that to develop the theory of the matrix 
algebraic equations, quadratic or of higher order, might provide useful tools 
to overcome technical mathematical obstacles encountered in the resolution 
of physical problems both in a classical and in a quantum context. 
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A APPENDIX 

Let's give few examples clearly illustrating our affirmation that the fun- 
damental theorem of algebra does not hold for matrix equation; thus, our 
examples point out QMEs having no solvent or infinitely many solvents as 
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well as equations with one or finitely-many solvents. In fact, these differ- 
ent scenarios are mainly traceable back to the existence of the square root 
of a n X n matrix, say M, which is a multivalued operation generally not 
representable in the form of a polynomial in M (Horn 1999). 

Firstly, let's consider the quadratic matrix equation X'^ — 2EX — Cq = 0, 
where Cq is the square matrix of order n with its diagonal entries equal to 
— 1, the elements of the first superdiagonal equal to unit and all the other 
vanishing. This equation may be immediately written as (X — E)'^ = Cq + E 
and it does not admit solvents since the square roots oi Cq + E = H^"^' does 
not exist (Gantmacher 1998). 

The matrix equation X^ -|- EX = Cq with Cq n x n diagonal matrix 
having (i2o)ii = 2(5i,i admits finitely-many diagonal solvents in the form 



5*1 = {1, €2, £3, . . . , e„}, S-2 = {-2, 1/2, 1^3, • • • , J^n} 



(43) 



where each of the diagonal entries €2, £3, . . . , e„ in Si and z^2, ^'3, • • • , i^n in 
S-2 may arbitrarily assume the values or —1. 

Finally, the matrix equation X^—EX = 2E provides an example wherein 
infinitely-many solvents exist having the form of lower triangular matrices 



X 



X 



( -1 
ci 2 


C2 



Sn/2 





ci 


C2 



Cn-1/2 





2/ 



n—even 



(44) 



-1 / 



n—odd 



where all the non-indicated entries are equal to and ci, C2, . . . , cn (cn^) G 

2 2 
C 
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B Appendix 

The n-quadratic matrices 



Si 



( 1 









2 \ 

2 

1 / 



, 5*2 



/ -1 
-2 -1 

\ -2 -2 -2 




-1/ 



(45) 



cannot be both solvents of the same QME due to fact that the matrix 
equation which furnishes the coefficient £i 



£1(5*1 — 5*2) — {Si — S2) 



(46) 



admits no solution. The incompatibility of equation is based on the Rouche- 
Capelli's theorem because rank{Si — S2) / rank{Sf — 5*1), where 



5*1 — <S'2 



/ 



sf - Si 




-1-22 



1-22 




V -(n-l)-22 -(n-2)-22 
Further, the n- quadratic matrices 



/ 2 

2 

2-22 
1-22 



(n - 3) • 22 



2 \ 
2 



(47) 



2 ... 2 / 

. {n-l)-2^ \ 
. (n - 2) • 2^ 





S[ 



( 1 

1 



V 



\ 



1 
3/ 



S'o 



)'^2 



/ 1 

1 



V 



are two solvents of the class of the QME for which 



/ an 
021 



Ci 



ai2 

0.22 



Oln-l 
0'2n-l 



\ 



1 
2/ 



\ 




(48) 



\ a„i a„2 
/ 1 — an — «12 ••• — «in-i 

«21 1 — 022 ••• — «2n-l 



«?i-ln-l 



^0 



-On-ll 



-an-12 
— «n2 



J- On—ln—l 
Otnn—l 



O-nn— 1 

\ 




-6/ 



5 / 



, (lij t 



(49) 
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